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Method and System for Automated Location-dependent Recognition of 

Flood Risks 

This invention relates to a method and system for automated 
location-dependent recognition of flood risks, flood states being measured and 
s location-dependent probability values being determined. In particular, the 

invention relates to a system and a method where a country-specific flood zone 
table, such as e.g. the First American 1 00-year flood zone table, is used to 
derive high resolution data on vulnerability factors for generalized insurance risk 
factors. 

i0 Today an appropriate automated flood rating system is painfully 

lacking. For many countries, it is hardly possible to do a technically correct flood 
rating and/or determination. A glance at the loss history shows that economic 
losses caused by flood are eaqually high or higher than those of earthquakes, 
windstorms other perils. For most of these other perils mostly exist already 

is various rating and/or prediction and/or early warning systems. High amounts of 
money, industrial power and time are lost by flood events. Additionally, with the 
trend of increasing insurance penetration for floods, the insurance and re-> 
insurance industry is affected ever more by flood losses. To extend the early 
warning and flood rating to detailed and even facultative business, however, the 

20 threat of immense data amounts has to be coped with. This is done by 
completely new concepts on the hazard event set as well as on the system 
and/or method side. Figure 2 and 3 show two diagrams which illustrates in the 
observation period 1 January 1997 - 31 December 2000; flood accounts for 
55% of all natural peril-induced damage to property cooperate business (Figure 

25 2) and for 44% of all natural peril-induced damage to industry business (Figure 
3). The reference numeral A/E shows flood related losses, B/F storm related 
losses, C/G earthquake/volcano related losses and D/H other losses. 

One of the objects of the invention is to provide a new approach to 
creating a stochastic event set. The main benefit of the new approach has to be 
30 the economic use of storage, which also positively influences the performance 
of the rating. Instead of storing water depths for each of the tens of thousands 
of events, e.g. only a flood frequency map or something similar should be 
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stored at a highly detailed level (90m resolution). The event return periods 
should be generable and storable on an aggregated high-resolution grid level 
for the territory. 

It is well known that floods cause equally high or higher economic 

5 losses than storm or earthquakes. Considering losses on industrial books, flood 
losses typically exceed by far earthquake or storm losses. Figure 1 shows in the 
observation period 1 Januaryl 997 - 3 1 December 2000, that flood accounts for 
55% of all natural peril-induced damage to property corporate business and for 
44% of all natural peril-induced damage to industry property business of a 

10 specific reinsurance company. This concerns not only industry itself, but with 
the permanent rise of the flood insurance penetration, in addition to a 
concentration of values in areas prone to flood, the insurance and re-insurance 
industry is affected ever more by flood losses. However, flood rating has been 
badly neglected, despite the fact that often detailed information (address) on 

15 the ceded risks is available. This is mainly owing to the fact that no efficient 
flood rating system is available. Ideally, such a system is made up of two 
components: first, it needs a street address geo-encoding to translate the 
addresses into a lat/lon and to make a lookup on the existing flood zones. 
Second, based on the lat/lon and flood zone information, a fully technical risk 

20 rating (including event-based capacity calculation) can be performed. 

The advantages of being able to do a correct technical risk rating for 
flood (e.g. on an Aggregated Loss Modelling (ALM) or Detailed Loss Modelling 
(DLM) basis) are plenty: (i) Better risk quantification: both components of an 
automated flood risk system will allow better estimatation of the risk exposure of 

25 a detailed property schedule. It will be possible to reduce the number of 
"unexpected losses" e.g. to insurance, reinsurance or other industry. The 
improved risk quantification is enhanced by taking into account the TC-Flood 
correlation; (ii) Better structuring of flood-exposed programs or portfolios: since 
the contribution of every single risk to the expected losses can be quantified, it 

30 is possible to determine site-dependent sub-limits or even to exclude certain 
sites from flood coverage; (iii) Essential improvements in flood risk assessment: 
It will be the first time in industry that a detailed and fully stochastic flood rating 
system is achieved. This can have a positive effect on the positioning of 
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industry and/or insurance and re-insurance in the market, and can generate 
new business opportunities; (iv) The availability of a detailed automated flood 
rating system will allow better qualification of property schedule data, e.g. in 
relation to identifying which policies have a flood coverage with the insurance 
5 industry; (v) Improved determination of capacity requirements: With a new 

event set approach it will be possible to allocate expected losses consistently to 
a specific contract. 

In the state of the art the US Geological Survey (USGS) National 
Elevation Dataset (NED) has been developed by merging the highest- 
. 10 resolution, best-quality elevation data available across the United States into a 
seamless raster format. NED is the result of the maturation of the USGS effort 
to provide 1:24,000-scale Digital Elevation Model (DEM) data for the 
conterminous US Projection. NED has a resolution of one arc-second 
(approximately 30 meters) for the conterminous United States. In the NED 

is assembly process the elevation values are converted to decimal meters as a 
consistent unit of measure, North American Datum of 1983 (NAD83) is 
consistently used as horizontal datum, and all the data are recast in a 
geographic projection. NAD83 is based on the Geodetic Reference System of 
1980 (GRS 80) ellipsoid that was developed to best-fit sea level over the entire 

20 earth. In the state of the art there can also be found the National Landcover 
Dataset (NLCD). The National Land Cover Characterization project was created 
in 1995 to support, the original Multi-Resolution Land Characterization (MRLC) 
initiative and fulfill the requirement to develop a nationally consistent land cover 
data set from MRLC data called National Land Cover Data 1992 (NLCD 92). 

25 This culminated in the September 2000 completion of land-cover mapping 
using a modified Anderson level II classification for the conterminous United 
States. The source data is the Landsat TM data for the years close to 1992. 
The National Land Cover Data (NLCD) is a2l-class land cover classification 
scheme applied consistently over the United States. In addition to satellite data, 

30 scientists used a variety of supporting information including topography, 

census, agricultural statistics, soil characteristics, other land cover maps, and 
wetlands data to determine and label the land cover type at 30 meter 
resolution. The spatial resolution of the data is also in geographic projection, 
NAD 83. The relevant classes for flood assessment are the following: A) Water: 
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(11) Open Water, (12) Perennial Ice/Snow; B) Developed: (21) Low Intensity 
Residential (Discontinuous Urban) including areas with a mixture of constructed 
materials and vegetation, where constructed materials account for 30-80% of 
the cover and vegetation for 20-70%. These areas most commonly include 

5 single-family housing units and population densities lower than in high intensity 
residential areas. (22) High Intensity Residential (Continuous Urban) including 
highly developed areas where people reside in high numbers. Examples 
include apartment complexes and row houses. Vegetation accounts for less 
than 20% of the cover and constructed materials for 80-100%. (23) 

io Commercial/lndustrial/Transportation (Industrial) including infrastructure (e.g. 
roads, railroads, etc.) and all highly developed areas not classified as High 
Intensity Residential. The growing need for land cover and other geospatial 
data within the federal government culminated in reforming the MRLC 
Consortium again in 2000 (MRLC 2001 ), with a second generation national 

15 land-cover data set to be developed from the data called National Land Cover 
Data 2001 (NLCD 2001). The production of NLCD 2001 will be implemented in 
a phased approach using mapping regions developed by USGS. Full 
production development is anticipated to begin in 2002, with completion 
targeted for 2004. Unfortunately, this data will not be available before 2005. 

20 Further there is the National Hydrography Dataset (NHD) in the state 

of the art, which is a set of digital spatial data that contains information about 
surface water features such as lakes, ponds, streams, rivers, springs and wells. 
The National Hydrography Dataset (NHD) is a newly combined dataset that 
provides hydrographic data for the United States on a hydrologic unit basis (8- 

25 digit number) . It is currently based on the content of the USGS 1:100,000-scale 
data, giving it accuracy consistent with those data. The horizontal positional 
accuracy is such that 90 percent of points are within 50 meters of their true 
position. The locations of points, lines, and the boundaries of areas are 
encoded using geographic (longitude-latitude) coordinates. The horizontal 

30 datum is the North American Datum of 1983 (NAD83). In the state of the art the 
FloodMap from First American provides the most accurate, up-to-date flood 
hazard zones for risk assessment, insurance and re-insurance. This digital map 
database is based on Flood Rate Insurance Maps (FIRMs) prepared by the 
Federal Emergency Management Agency to depict flood risks in support of the 
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congressionally mandated National Flood Insurance Program (NFIP). The 
result of engineering flood insurance studies, FIRMs are used to implement 
floodplain management measures to reduce flood risks to communities and to 
perform flood insurance policy rating. These FIRMs depict Special Flood 

5 Hazard Areas that are subject to inundation for both the 100-year flood 

(national standard for flood insurance) and the 500-year flood zones. The First 
American FloodMapSM database has been cleaned and corrected to reflect 
information as shown on these paper FIRMs. National, State, Territory and 
County First American Flood Hazard Certification currently has 3*100+ counties 

10 of digital Flood Insurance Rate Maps. The county level coverage represents 
over 99% of the U.S. population. FloodMap reflects all FEMA revisions (over 
12000 a year). The entire database is refreshed quarterly to ensure accuracy, 
currency and investment value. The dataset is provided in geographic 
coordinates based on NAD83. In the dataset the coastal floodways and 

i5 floodplains (flood zone code "V*")are considered in a separate storm surge risk 
evaluation. Finally, there is a data source for the USGS Daily Values 
publication, the so-called USGS WATSTORE database in the state of the art. It 
includes stream flow and related stream data with remarks (i.e. stream flow, 
stage, water temp., lake level, reservoir storage, level below surface and more). 

20 It includes 38,972 stations covering a record period of maximally from 1781 to 
the present. It contains raw data, daily, monthly and exceedance statistics. The 
datum for the coordinates of the gauging stations (6) is also NAD83. 

It should be pointed out that, besides the method according to the 
invention, the present invention also relates to a system and a computer 
25 program product for carrying out this method. 

In particular the objects are achieved through the invention in that for 
automated location-dependent recognition of flood risks, flood states being 
transmitted to a central unit and location-dependent flood probability values 
being determined, in that the central unit comprises a multi-dimensional lookup 
30 table corresponding to a spatial high resolution grid based on decentralized 
measurements of flood risk factors of a specific territory whereas the flood risk 
factors (P) being associated with the grid indicating the averaged flood 
frequency and/or susceptibility to flooding within a grid cell, in that the system 
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comprises distributed gauging stations, whereas river discharge parameters (T) 
are measurable by the distributed gauging stations within a grid cell and 
transmitted over a network to the central unit, the river discharge parameters 
comprising at least values for return period and/or intensity of measured events, 

5 in that the central unit comprises a correlation-module generating an event- 
specific averaged probabilistic water depth value (H) for an flood event based 
on the linked flood risk factors and the river discharge values, and associating 
the probabilistic water depth value (H) to the corresponding grid cell. The 
system can comprise an cell arbitrator module accting on at least on grid-based 

10 composition module according to the avaraged probabilistic water depth values 
(H). The grid-based composition module can e.g. comprise at least an early 
warning system signaling flood risk in the appropriate cell and/or an automated 
damage prediction systems and/or damage covering systems for land-based 
installations based upon at least the flood risk factors. The grid-based 

15 composition module also can comprise at least an optimization and/or control 
module for protection installation of technical and/or industrial facilities based 
upon at least the flood risk factors. 

Under spatial high resolution a resolution with a cell size below 1km 2 
is understood, in particular a cell size below 10'OOOm 2 . One of the advantages 

20 of the invention is that a much more economic use of storage can be achieved, 
which also essentially influences the performance of the rating. For example, 
instead of storing water depths for each of tens of thousands of events or more, 
only one flood frequency map is stored at a highly detailed level (90m 
resolution). The event return periods are generated and stored on cell level. 

25 Both pieces of information are combined to calculate the water depth at a 
location for all the events. Additionally, a better structuring of flood-exposed 
programs or portfolios can be achieved. Since the contribution of every single 
risk to the expected losses can be quantified, it is possible to determine site 
dependent sub-limits or even to exclude certain sites from flood coverage. 

30 Further essential improvements in flood risk assessment can be provided: It will 
be the first time in industry that a detailed and fully stochastic flood rating 
system is achieved. This can have a positive effect on the positioning of 
industry and/or insurance and re-insurance in the market, and can generate 
new business opportunities. Another advantage of the invention is the 

35 availability of a detailed automated flood rating system, which will make it 
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possible to require better quality property schedule data, e.g. in relation to 
identifying which policies have a flood coverage with the insurance industry. 
Finally, improved capacity allocation and determination can be done: With a 
new event set approach, it will be possible to allocate expected losses 
5 consistently to a specific contract. The correlation module can e.g. comprise at 
least five adaptable correlation parameters Xi, x2,..., xs, whereas the adaptable 
parameters can e.g. be correlated by 

H = max 2 3 — ^ U 1 , 

For flood cell zones along a coast located at different heights in 
10 meters above sea level, different correlation modules can e.g. be used. For 
example, for flood cell zones along a coast located lower than 10 m above sea 
level can be determined by a specific correlation module. The determination of 
the flood cells along a coast can e.g. be based additionally upon storm surge 
events. The determination of the flood cells along a coast can also comprise 
15 e.g. a Sea, Lake and Overland Surges from Hurricanes (SLOSH) method. 

In an embodiment of the invention, a linking module comprises at 
least one adaptable event factor providing the spatial and/or temporal 
correlations for the discharge measurements of different gauging stations 6. 
This embodiment has inter alia the advantage that the distribution of return 
20 periods between the single cells can show spatial and temporal correlations 
determined from discharge measurements of gauging stations 6, which can be 
taken into account by the at least one adaptable event factor. 

In another embodiment of the invention, the river discharge 
parameters (T) are determined by measurements in a low spatial resolution in 

25 relation to the grid of the flood risk factors. This embodiment has inter alia the 
advantage that the use of storage can further be economized, which also 
essentially influences the performance of the rating. For instance, instead of 
storing water depths for each of the tens of thousands of events or more, only 
one flood frequency map is stored at a highly detailed level (90m resolution). 

30 The event return periods are generated and stored on cell level. Both pieces of 
information are combined to calculate the water depth at a location for all the 
events. 
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In still another embodiment according to the invention, the resolution 
of the grid cells is below 10boom 2 for the diameter. It must be pointed out that, 
in principle, it is possible to set the resolution to an arbitrary small cell size. 
However, the resolution is limited e.g. by the inherent errors of the parameters. 
5 This embodiment has, among other things, the advantage that the desired flood 
parameters can be properly determined by the system for a specific cell and/or 
site size. 

In a further embodiment according to the invention, the flood risk 
factors (P) are derived based upon at least geomorphologic parameters. The 

10 geomorphologic parameters comprise e.g. horizontal distance and/or elevation 
difference to the next river and/or drainage area. This embodiment has inter 
alia the advantage that it comprises a method which is known in the state of the 
art, and e.g. errors can be easily estimated. So the flood risk factors (P) can 
e.g. be derived by an interpolation module based upon a country-specific flood 

15 zone table depending on horizontal distance and/or elevation difference. For 
the United States of America the country-specific flood zone table can 
comprise, for example, the First American 100-year flood zone table. 

In an embodiment according to the invention, the vulnerability factors 
are determined based upon historical dataset of corresponding portfolios, and a 
20 generalized insurance risk is automated derivable from the vulnerability factors. 
This embodiment has, among other things, the advantage that a fully 
automated system for determining the vulnerability factors of places and sites 
can be provided. This has not been possible up to now in the state of the art. 

In another embodiment according to the invention, the system is 
25 accessable over a network by client nodes, whereas the system comprises a 
billing module with a billing gateway interface for access to the central unit first 
call datail records of a client node being transmittable from the central unit to 
the billing module. By means of a proxy module of the system second call detail 
records of the client node also can be downloadable from the central unit by 
30 means of the proxy module at least the identity of the client node and/or 
duration of the access to the central unit and/or service being able to be 
captured and able to be passed on to the billing module. By means of the billing 
module of the system TAP files corresponding to the obtained service are able 
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to be generated, an these are transmittable, together with billing instructions, to 
a clearing module, the billing instructions including at least user-specific and/or 
service-provider-specific billing data. This embodiment has, among other 
things, the advantage that a full billing of the use of the system for third parties 
5 can be provided. 

Additional features and advantages will become apparent to those 
skilled in the art upon consideration of the following detailed description of 
illustrative embodiments exemplifying the best mode of carrying out the method 
as presently perceived. 

10 The present disclosure will be described hereafter with reference to 

the attached drawings which are given as non-limiting examples only, in which: 

Figure 1 gives an shematic overview fo an embodiment of the 

invention. 

Figure 2 and 3 show two diagrams which illustrates in the 
15 observation period 1 January 1997 - 31 December 2000; flood accounts for 
55% of all natural peril-induced damage to property cooperate business (Figure 

2) and for 44% of all natural peril-induced damage to industry business (Figure 

3) . The reference numeral A/E shows flood related losses, B/F storm related 
losses, C/G earthquake/volcano related losses and D/H other losses. 

20 Figure 4 shows a diagram illustrating a flood frequency map (p-map) 

based on extended FEMA 100-y flood zone. The interpolation depending on 
horizontal and vertical distance to next river to get other return periods. The 
reference numeral 1 is a 50-y flood zone, 2 is a 100-y flood zone, 3 is a 500-y 
flood zone and 4 is a river. 

25 Figure 5 shows a diagram illustrating an example of flood frequency 

map, where as the spatial resolution is 90 m grid cell size for the chosen map. 

Figure 6 is a diagram which illustrates gauging stations 5 with 25 
years of records used for event set generation. This example shows used 
gauging stations 5 in North Caroline (US) where the determination is done 
30 based upon daily stream flow records. 
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Figure 7 is a diagram which illustrates the example event set with 
return periods per cell for the biggest event in North Carolina. 

Figure 8 shows a diagram which illustrates the relation between flood 
frequency (P-value), water depth (H) and even set (return period T). The 
5 relation is plotted on Fig. 10 with the black dots indicating the values related to 
100-year and 500-year flood zones and 100, respectively 500-year flood 
events. 

Figure 9 shows a diagram which illustrates the same relationship as 
in Figure 8 but with different axes, representing the relationship between flood 
10 exposure, return period and water depth. It can be noted that for this example a 
flood event with a return period of 10 years, can only affect locations with P- 
values higher than 0.82 with water depths ranging up to 0.8 m. For a 100-year 
flood this range varies from 0.0 m (P<0. 54) up to 2.2 m (P=1 ). 

Figure 10 shows a diagram illustrating the loss in relation to the 
15 frequency of the event for a random portfolio 6 of industry, residential and/or 
buildings etc. and a wind portfolio on the example of North Carolina (with 
vulnerabilities). 

Figure 1 1 shows an elevation model of North Carolina (US) with the 
zones mentioned below. The elevation is given in meter above sea level (a.s.l.). 

20 The overall operation, method and system of one embodiment of the 

subject method is discussed in the following. The reference numeral 20 in 
Figure 1 refers to a central unit 20 having the necessary infrastructure, 
including hardware and software components and/or units, to achieve a 
described system and/or method according to the invention. The system and 

25 method allow creation of a stochastic event set. One of the main benefit of the 
new approach is the highly economic use of storage, which directly influences 
the performance of the system and its operation. Instead of storing water 
depths for each of tens of thousands of events, only one or very few flood 
frequency maps have to be stored at a highly detailed level (90m resolution). 

30 The event return periods are generated and stored on the desired cell level. 
Both pieces of information are combined to calculate the water depth at a 
location for all the events. In this embodiment variant, a central unit 20 
comprises a multi-dimensional lookup table 203 corresponding to a spatial high 
resolution grid 60/61 based on decentralized measurements of flood risk factors 
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of a specific territory whereas the flood risk factors (P) being associated with 
the grid indicating the averaged flood frequency and/or susceptibility to flooding 
within a grid cell 60/61 . As an example of the flood risk factors, i.e. a flood 
frequency map, e.g. the state of North Carolina (NC) / United States (US) can 

5 be chosen as prototype territory. North Carolina has a surface of 140'0OO km 2 
(for comparison mainland United Kingdom has 220 000 km 2 ) and 6.8 Mio 
citizens. The extent of North Carolina is 850 km in west-east direction and 305 
km in north-south direction. The highest point culminates on Mount Mitchell at 
2037 meters above sea level. The Atlantic coast is 484 km long. The 

10 topography is divided into three distinct regions: the Atlantic coast plain (2/3) up 
to 80 km into the inner land (elevation < 90 m a.s.l.), the Piedmont Plateau 
(elevation 90-460ma.s.l.) and the Appalachian (elevation > 460 meters above 
sea level). Figure 11 shows an elevation model of NC with the mentioned 
zones. Between Cape Hatteras and Cape Fear some regions are often affected 

15 by hurricanes. For the flood frequency map (or P-map) the grid can e.g. consist 
of values between 0 (never flooded) and 1 (always flooded) indicating the 
susceptibility to flooding. In this example the values can e.g. be derived from 
the First American 100-year flood zones, the digital elevation model and the 
river network. As mentioned, the flood frequency map, i.e. the flood risk values 

20 associated with the grid, can e.g. consist of values between 0 (never flooded) 
and 1 (always flooded) indicating the susceptibility to flooding of each location 
(latitude, longitude). The values of the flood risk factors (P) between 0 and 1 
can be related directly to flood frequency, for example P=0.58 stands for the 
100-year flood zone. Of course, other flood zones could be chosen as e.g. flood 

25 zones with a certain other return period (50-, 100-, 250-, 500-years). The flood 
zones can be determined based on geomorphologic parameters such as the 
horizontal distance and elevation difference to the next river and/or the 
drainage area. The data in this example were e.g. validated with FEMA data 
from the USA and proofed good results in several countries (CZE, UK, SVK, 

30 BEL). Figure 4 shows a diagram illustrating a flood frequency map (p-map) 
based on extended FEMA 100-y flood zone. The interpolation depending on 
horizontal and vertical distance to next river to get other return periods. The 
reference numeral 1 is a 50-y flood zone, 2 is a 100-y flood zone, 3 is a 500-y 
flood zone and 4 is a river. In the example of the embodiment, the First 

35 American 100-year flood zones are used to interpolate by an interpolation 
module depending on the horizontal distance and elevation difference of the 
cell. The interpolation module can be achieved through hardware and/or 
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software. An additional advantage of taking the 100-year flood zones for 
calibration is that flood protection measures will be automatically taken into 
account. The processing steps can comprise e.g. the following steps: (i) 
Prepare National Elevation Dataset: Merge, project and fill grid; (ii) Prepare 

5 First American FloodMapSM Dataset: Select zones, project, clean and expand 
grid by 600 m (the size of the buffer will be adapted depending on the regions), 
which results in a extended or buffered 100-year flood zone, (iii) Prepare 
National Hydrography Dataset: Select rivers, project grid. (All the datasets 
contain the elevation data and can e.g. be saved to a storage module); (iv) 

10 Generation of p-value map: With a appropriate processing module the input 
files are then processed to result in a grid with the extend of the buffered flood 
zones with p-values. a) Nearest neighbor lookup of buffer grid points on river 
network: Determination of horizontal distance and elevation difference, b) 
Nearest neighbor lookup of grid point at 100-year flood zone border on river 

15 network: Determination of horizontal distance and elevation difference, c) 
Determine calibration parameter with MARS formula, d) Smoothing of 
calibration parameter, e) determine p-value and e.g. save them into to storage 
module, (v) Post-processing: Aggregate, clean and reproject grid to geographic 
- coordinates. 

20 The quality of the p-value map can e.g. be checked by a test 

portfolio with known locations. These locations can then be checked with both 
grids (First American flood zone and p-value map). The results for this example 
are illustrated in Table 1, where all locations lying within the First American 100- 
year flood zones are taken and checked against the p-value map. As can be 

25 seen, 73% of all checked locations in North Caroline (NC) are also found in the 
100-year flood zone of the p-value map, and 97% are in 500-year flood zones 
of the p-value map. It is important to note that the flood zone along the coast 
can often include storm surge zones, and it is clear that in these areas the 
presented method is not well adapted (no rivers available). Therefore a second 
f 30 check can be made, taking only the locations above 10 meters above sea level, 
and in that case 92% of the locations from the p-value map fall into the 100- 
year flood zone. 

Table 1. Comparison of locations lying in First American 100-year 
flood zone with those from p-value map. 
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Flood zone 


All locations within First 
American 100-vear flood 
zone 


Locations within First 
American 100-vear flood 
zone (>10 m. a.s.l.) 


50-year 


73 (47%) 


39 (74%) 


100-year 


112(73%) 


49 (92%) 


250-year 


146 (95%) 


52 (98%) 


500-year 


149 (97%) 


53(100%) 



Over the whole of NC, the area affected by a 100-year flood 
determined by the system from the p-value map underestimates the First 
American flooded area by only 2.6%. Comparing the two datasets, it can be 
seen for this example that in the flat coastal region the flooded area is rather 
underestimated, which is due to the fact that storm surge is not taken into 
account yet. A slight overestimation exists in the central part of NC, but this is 
not significant. The final p-value map has in this embodiment example a cell 
size of 90 m, which is precise enough taking into account the precision of the 
geo-encoding. By that way, the total amount of data can be divided by a factor 
9. The total size of the p-value map (p>0.2) for NC includes 5.2 mn cells. 
Taking the whole USA, this would result in 360 Mio cells, assuming that the 
same percentage of area is flooded. 

Further the system comprises distributed gauging stations 
6/30/31/32. River discharge parameters (T) are determined by the distributed 
gauging stations 6/30/31/32 within a grid cell 60/61 and transmitted over a 
network 50 to the central unit 20, the river discharge parameters comprising at 
least values for return period and/or intensity of measured events. For this 
example the river discharge parameters (T) are also called event set T-table. 
The T-table comprises for the example the return periods of 1183 flood events 
affecting the 1082 cell zones of North Carolina. A probabilistic event set is 
derived from the statistics of 74 river discharge gauging stations 6/30/31/32. 
Again, the event set specifies the return period or intensity of one event in a 
given cell of the grid. The distribution of return periods between the cells can 
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take into account the spatial and temporal correlations determined from 
discharge measurements of gauging stations 6/30/31/32 covering a range of 25 
years from 1977-2002. Totally, 74 stations fulfilled this requirement for the 
example of North Carolina. For the whole USA, approximately 3000 stations 
5 can be used. It is clear that not only head catchment gauging stations 

6/30/31/32 can be used for the system, but also stations at main rivers. It has to 
be noted that only few stations were available in the coastal area as can be 
seen on Figure 6. 

The reason for the selection of return period as a measure of 
10 intensity lies mainly in its relative weak variation along the river network, i. e. the 
return period between adjacent cells of the grid on a river varies "smoothly". 
Figure 7 shows the distribution of return periods on cell resolution for the 
biggest event in the North Carolina event set. The chosen definition of events 
through return periods further allows real events to be easily included by the 
15 means of extreme value analysis of the time series, for example the 1993 flood. 
Real events also provide the possibility, on the one hand, of a direct 
comparison with possible loss data and, on the other hand, a plausibility check 
of the probabilistic events. Finally, post-event analysis could be made by 
downloading the newest discharge measurements and analyzing them with the 
20 past time series. 

The central unit 20 further comprises a correlation-module 2 1 . An 
event-specific averaged probabilistic water depth value (H) for an flood event 
based on the linked flood risk factors and the river discharge values is 
generated by the correlation module 2 1 , and associating the probabilistic water 

25 depth value (H) to the corresponding grid cell 60/61 . The event-specific 

averaged probabilistic water depth value (H) for an event is in this example the 
desired location-dependent flood probability value. The correlation module 21 
can e.g. be achieved through hardware and/or software. To link the flood 
frequency map (p-value map) and the event set (return period T)with the 

30 hazard, i.e. water depth H, a relationship indicating a mean average water 

depth as a function of these two parameters is necessary. To determine such a 
relationship e.g. an analysis based on the probabilistic flood event set of a 
known territory can be done. In this example, a probabilistic flood event set of 
the UK was chosen. The chosen event set consists of 973 flood maps 
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indicating the water depth at each flooded cell with a resolution of 50 m. For the 
analysis, the unprotected event set was used. In order to perform the analysis, 
the flood frequency P-value of each cell had to be known. Therefore, the P-map 
determined by the system was taken as a mask. Specifically, this means that 
for each event or flood footprint, the water depths of all cells with a certain P- 
value were collected. This collection allows e.g. for a classification into 8 water 
depth classes (0-4.0 m) and 16 p-value classes 0.2-0.9. If the number of 
occurrences in each class is divided by the total sum of observations, the 
relative frequency f, can be calculated as 



This method can also be explained by the fact that, with a high 
number of events and flooded cells, this frequency distribution of water depths 
and P-values becomes an intrinsic function which is independent of the number 
of events. Or, in other words, if you add 10 additional events, the relative 
frequency will not change, and still represents the occurrence frequency. 

To give an analytical expression for the relationship, the values of 
the example were fitted to the following function 



f = 



_NumbeofoccurrencsBnclass/_ 
^ NumbeDfoccurrencsanclassT 



and the corresponding return period is 




T=expe.9+0.5H-8.0 P + 1.2 H P) 



or 




The function is plotted on Figure 8 with the black dots indicating the 
values related to 100-year and 500-year flood zones and 100, resp. 500-year 
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flood events. The P-values for the 100-year (P>0. 58) and 500-year flood zone 
(P>0.35) were given by the calibration of the MARS model to the FEMA flood 
zones of several regions in the USA. Taking the P-H-T relationship, these P- 
values with the corresponding return period T result in water depths of 0.3 m, 

5 resp. 0.1 m. In Figure 9 the same relationship is presented with different axis. It 
can be noted that, for example, a flood event with a return period of 10 years 
can only affect locations with P-values higher than 0.82, with water depths 
ranging up to 0.8 m. For a 100-year flood, this range varies from 0.0 m (P<0. 
54) up to 2.2 m (P=1 ). To get additional information for this example about the 

10 mentioned relationship, a questionnaire was used where all flood people from 
Zurich, Munich and Armonk gave their best estimate of such a relationship 
between flood frequency, event return period and water depth. The results of 
this "expert view" are compared with the analytical relationship. The comparison 
shows that the analytical relationship is slightly more conservative than the 

15 estimates indicating that during the calibration work some smaller modifications 
to the relationship could be done. 

It is clear that the system and method also allow the values to be 
verified and adapted for flood zones other than 100-year and 500-year. The 
proposed P-values can be e.g. as follows: (i) 50-year flood zone: P=0.69 

20 (resulting water depth with 50-year flood event H=0.4 m); (ii) 250-year flood 

zone: P=0.45 (resulting water depth with 250-year flood event H=0.2 m); (iii) 25- 
year food zone: P=0.82 (resulting in water depth with 25-year flood event H=0.6 
m). Based on the return periods of all events (in every second cell) and on the 
complete flood frequency map, the water depth distribution can be determined. 

25 In the example given here in 54% of all cases the water depth is below 50 cm. 

As shown above for the example, a rather large area of North 
Carolina lies below 10 meters above sea level, and is potentially affected by 
storm surge events. The system and method described above is well adapted 
for riverfloods inland located higher than 10 meters above sea level, but along 
30 the coast the flood zones also include the risk of storm surge, and for this area 
it can be reasonable to use a specific model. Such differentiations can easily be 
included in the same way, and form part of the invention. For example, storm 
surge events can be included in the system and/or method. Another possibility 
would be to use the SLOSH (Sea, Lake and Overland Surges from Hurricanes) 
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method, which is a method used by the National Hurricane Center (NHC) to 
estimate storm surge heights resulting from historical, hypothetical, or predicted 
hurricanes. The second model takes into account pressure, size, forward 
speed, track and winds, which could be obtained from TCart. Independently of 

5 the method and/or system used, some kind of storm surge frequency map 

should be developed in order to be consistent with the river flood event set. The 
storm surge event set would then also consist of return periods, where the 
events are "linked" to the corresponding TC events. For storm surge different P- 
H-T relationships as well as vulnerability curves must be defined. The relation 

10 between river floods and TCs is less evident. However also the correlations 
between the river floods and TCs can be considered. Here, it can be 
reasonable to perform a matching between the TC rainfall maps and the river 
flood return periods. 

Further the system can comprise an cell arbitrator module 22 accting 
on at least on grid-based composition module 23 according to the avaraged 
probabilistic water depth values H. The grid-based composition module 23 can 
comprise at least an early warning system signaling flood risk in the appropriate 
cell 60/61 . Also, the grid-based composition module 23 can comprise 
automated damage prediction systems and/or damage covering systems for 
land-based installations based upon at least the flood risk factors. The grid- 
based composition module 23 can comprise at least an optimization and/or 
control module for protection installation of technical and/or industrial facilities 
based upon at least the flood risk factors. To consider correlations caused by 
correlations the system comprises a linking module with at least one adaptable 
event factor providing the spatial and/or temporal correlations for the discharge 
measurements of different gauging stations 6/30/31/32. The gauging stations 
6/30/31/32 comprise at least a sensor and/or measuring fixture to measure river 
water level parameters, based upon which river water level parameters the river 
discharge parameters are determined. 

30 The system can be accessable over a network 5 1 by client nodes 

40/41/42, in particular by mobile client nodes. The reference numerals 50 and 
51 accordingly stand for the various heterogeneous networks, such as e.g. a 
wired LAN, i.e. a local fixed network, in particular also the PSTN (Public 
Switched Telephone Network), etc., a Bluetooth network, e.g. for installations in 
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covered-over localities, a mobile radio network with GSM and/or UMTS, etc. or 
a wireless LAN. Therefore it is possible to integrate the system f.e. as a part of 
information and/or early warning system for third parties. To be understood as 
then client nodes 40/41/42 can be inter alia all possible so-called Customer 

5 Premise Equipment (CPE) intended for use at various network locations and/or 
in different networks. Therefore the client nodes 40/41/42, illustrated 
schematically in Figure 1, can be a mobile radio telephone 40 or a laptop- 42, 
palmtop 41 or PDA- 4 1 computer (Personal Data Assistant) equipped to 
communicate over the network 5 1 , in particular a mobile radio network. The 

10 client nodes 40/41/42 possess one or more different physical network 

interfaces, which can also support a plurality of different network standards. 
The physical network interfaces of the client nodes 40/41/42 can comprise e.g. 
interfaces for Ethernet or for another wired LAN (Local Area Network), 
Bluetooth, GSM (Global System for Mobile Communication), GPRS 

15 (Generalized Packet Radio Service), USSD (Unstructured Supplementary 

Services Data), UMTS (Universal Mobile Telecommunications System) and/or 
WLAN (Wireless Local Area Network), etc. The client nodes 40/41/42 can 
include an identification module, preferably a chipcard fixed or removably 
connected to the client nodes 40/41/42, preferably a so-called SIM-card 

20 (Subscriber Identification Module). The identification module contains an 
identification data element, for example an International Mobile Subscriber 
Identity (IMSI). Furthermore, the client nodes 40/41/42 can include a data store, 
preferably located on the chipcard, which can be used as an account to store 
monetary amounts. The system can comprise at least one billing module with a 

25 billing gateway interface for access to the central unit 20 first call datail records 
of a client node being transmittable from the central unit 20 to the billing 
module. By means of a proxy module of the system second call detail records 
of the client node 40/41/42 are downloadable from the central unit 20 by means 
of the proxy module at least the identity of the client node and/or duration of the 

30 access to the central unit 20 and/or service being able to be captured and able 
to be passed onto the billing module. By means of the billing module of the 
system TAP files corresponding to the obtained service can be generated, an 
these can then be transmitted, together with billing instructions, to a clearing 
module, the billing instructions including at least user-specific and/or service- 

35 provider-specific billing data. 
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Although the present disclosure has been described with reference 
to particular means, materials and embodiments, one skilled in the art can 
easily ascertain from the foregoing description the essential characteristics of 
the present disclosure, while various changes and modifications may be made 
5 to adapt the various uses and characteristics without departing from the spirit 
and scope of the present invention as set forth in the following claims. 
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